clear all
% load ParEst_eh0_time0_Mex_girls
% Chute = [1.32, ParEst_eh0_time0_Mex_girls(7:10)];
load ParEst_eh0_time2_Mex_girls
Chute = ParEst_eh0_time2_Mex_girls;
HandleLo = @Principal_A2_n;
LoBo = [-1.25, 0.1, 0.1, 0.1, 0.1];
HiBo = [2.75, 10, 10, 10, 10];
                 options = saoptimset( ...
                    'PlotFcn',{@saplotbestf,@saplottemperature,@saplotf,@saplotbestx,@saplotx},...
                     'TolFun',3000,'ReannealInterval',4, ...
                     'InitialTemperature',10*10^2*[4 10 10 10 10], ...
                     'TemperatureFcn', @temperaturefast, 'StallIterLimit',1000);
[x,fval,exitflag,output]  = simulannealbnd(HandleLo, Chute, LoBo, HiBo, options)
ParEst_eh0_time2_Mex_girls = x;
save ParEst_eh0_time2_Mex_girls ParEst_eh0_time2_Mex_girls